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SUMMARY 


A gas-chromagraphic  technique  has  been  employed  to  measure  the  solubilities 
of  diethylamlne  and  ethyl  acetoacetate  in  water  at  low  concentrations.  Diethyl-* 
amine  solubility  was  found  to  Increase  with  decreasing  concentration  and 
Increasing  acidity.  Ethyl  acetoacetate  solubility  was  found  to  decrease  with 
decreasing  concentration. 

Mathematical  expressions  were  formulated  to  describe  the  observed  solubility 
behavior,  and  these  were  incorporated  into  the  EPAEC  scavenging  model  to 
predict  concentrations  of  these  materials  in  rain  for  comparison  with  experi- 
mental measurements.  Agreement  between  experiment  and  theory  was  found  to 
improve  over  previous  estimates,  which  were  based  upon  less  accurate  solubility 
data.  Diethylamlne  results  showed  generally  good  agreement,  while  those  for 
ethyl  acetoacetate  exhibited  considerable  deviation.  Disagreement  in  this 
latter  case  was  attributed  primarily  to  the  neglect  of  aqueous-phase  mixing 
effects  in  the  model  calculations. 
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INTRODUCTION 

The  objective  of  the  research  described  In  this  report  Is  the  measurement 
of  the  water  solubility  of  dlethylamlne  and  ethyl  acetoacetate  at  low  con- 
centrations. Interest  in  this  Information  has  arisen  from  previous  Army 
research  (LEE  and  HALES  (1974))  dealing  with  precipitation  scavenging  of 
organic  compounds,  which  demonstrated  a marked  deviation  between  observed 
and  predicted  scavenging  rates. 

The  basis  for  predicting  scavenging  rates  in  this  former  work  was  the  EPAEC 
scavenging  code,  which  had  been  demonstrated  previously  to  predict  sulfur 
dioxide  scavenging  with  reasonable  accuracy  (DANA,  HALES,  and  WOLF  (1973)). 

An  analysis  of  the  results  Indicates  that  the  most  probable  reason  for  the 
large  disagreement  exhibited  in  the  caae  of  the  organic  materials  was  poor 
estimates  of  the  water  solubility  of  these  compounds,  which  were  required 
as  model  Input. 

If  solubility  Is  expressed  In  terms  of  the  equation 

P - H'c£  , (1) 

where  p (atm.)  Is  the  partial  pressure  of  vapor  in  the  gas  phase  associated 
with  an  equilibrium  aqueous-phase  concentration  (moles/ml),  then  the 
solubility  parameters  H'  (atm.  ml/mole)  can  be  predicted  as  shown  In 
Table  1 (LEE  and  HALES  (1974)).  Observed  washout  concentrations  are  com- 
pared with  those  predicted  using  these  solubility  estimates  In  Tables  2-6. 

Additional  calculations  have  been  performed  showing  that  much  better  agree- 
ment between  experiment  and  theory  can  be  achieved  using  higher  solubilities 
In  the  case  of  dlethylamlne,  and  lover  ones  In  the  case  of  ethyl  acetoacetate 
This  rather  weak  evidence  that  the  solubility  data  in  Table  1 Is  Inappropri- 
ate, Is  supported  by  the  highly  approximate  nature  of  the  calculations 
leading  to  these  values.  Estimates  for  both  compounds  were  based  upon 
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TABLE  1.  ESTIMATED  SOLUBILITIES  OF  ETHYL  ACETOACETATE 
AND  DIETHYLAMINE  IN  WATER  AT  5*C 


« 


Compound  H*  (atm,  ml /mole) 

Ethyl  acetoacetate  0.25 

Diethy  lamina  11. 


TABLE  2.  OBSERVED  VERSUS  PREDUCTED  WASHOUT  CONCENTRATIONS  - EEA 
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Rain 

Measured 

Predicted 

Rain 

Measured 

Predicted 

Collector 

Concentration 

Concenratlon 

Collector 

Concentration 

Concentration 

Station 

^obs. 

CH'-0.25 

Station 

C . 
oba. 

CH'-0.25 

A 3 

<0.10  ppm 

0.04  ppm 

B 8 

0.05  ppm 

4 

0.16 

0.18 

9 

— 

0.21 

5 

0.50 

0.71 

10 

— 

0.67 

6 

1.1 

2.41 

11 

1.0  ppm 

1.68 

7 

1.5 

7.96 

12 

0.70 

3.18 

8 

1.6 

27.7 

13 

1.2 

4.43 

9 

1.5 

13.9 

14 

0.78 

4.21 

10 

1.2 

4.48 

15 

0.73 

2.54 

11 

1.0 

1.12 

16 

0.33 

0.97 

12 

0.32 

0.25 

17 

0.35 

0.27 

13 

0.10 

0.05 

18 

0.12 

0.06 

TABLE  3.  OBSERVED  VERSUS  PREDICTED  WASHOUT  CONCENTRATIONS 

- EEA 

Rain 

Measured 

Predicted 

Rain 

Measured 

Predicted 

Collector 

Concentration 

Concentration 

Collector 

Concentration 

Concentration 

Station 

C . 
obs. 

Si’  — .25 

Station 

C . 
obs. 

CH'-0.25 

A 1 

<0.10  ppm 

0.27  ppm 

B 5 

<0.10  ppm 



2 

<0.10 

0.89 

6 

0.26 

3 

0.80 

2.38 

7 

0.16 

0.16  ppm 

4 

1.0 

5.23 

8 

0.31 

0.49 

5 

1.5 

9.99 

9 

0.36 

1.23 

6 

2.0 

28.7 

10 

0.31 

2.44 

7 

1.6 

9.99 

11 

0.27 

3.73 

8 

0.77 

5.23 

12 

0.37 

4.31 

9 

0.20 

2.38 

13 

0.28 

3.73 

10 

0.11 

0.89 

14 

0.12 

2.44 

11 

<0.10 

0.27 

15 

<0.10 

1.23 

12 

<0.10 

— 

16 

<0.10 

0.49 

13 

<0.10 

— 

17 

<0.10 

0.16 

18 

<0.10 

TABLE  4.  OBSERVED  VERSUS  PREDICTED  WASHOUT  CONCENTRATIONS  - DEA 


ft 


ft 


ft 


ft 


ft 


ft 


ft 


ft 


Rain 

Collector 

Station 

Measured  Predicted  Rain  Measured 

Concentration  Concentration  Collector  Concentration 

C , Cu.  . . Station  C . 

obs.  H'«  11  obs. 

Predicted 

Concentration 

V-  11 

A 4 

0.6  ppm 

— 

B 11 

0.9  ppm 

0.002  ppm 

5 

0.7 

— 

12 

1.7 

0.021 

6 

2.1 

0.005  ppm 

13 

2.2 

0.112 

7 

1.4 

0.042 

14 

2.7 

0.304 

8 

1.3 

0.225 

15 

3.0 

0.425 

9 

2.1 

0.628 

16 

2.1 

0.304 

10 

2.5 

0.898 

17 

2.8 

0.112 

11 

1.0 

0.628 

18 

1.1 

0.021 

12 

0.6 

0.225 

19 

1.2 

0.002 

13 

0.1 

0.042 

14 

0.1 

0.005 

TABLE  5.  OBSERVED  VERSUS  PREDICTED  WASHOUT  CONCENTRATIONS  - DEA 


Rain 

Collector 

Station 

Measured  Predicted  Rain 

Concentration  Concentration  Collector 

Cob..  °H'.  11  St*tlon 

Measured  Predicted 

Concentration  Concentration 

Cobs.  Si'  ■ 11 

A 6 

0.19  ppm 

— 

A 12 

0.21  ppm 

— 

7 

0.55 

13 

0.42 

— 

8 

0.62 

0.014  ppm 

14 

1.3 

— 

9 

0.96 

0.357 

15 

1.3 

0.013  ppm 

10 

1.3 

1.4 

16 

1.7 

0.19 

11 

1.0 

0.357 

17 

2.2 

0.48 

12 

0.85 

0.014 

18 

1.7 

0.19 

13 

0.56 

— 

19 

1.4 

0.013 

14 

0.39 

— 

20 

1.8 

21 

1.3 

22 

0.80 

23 

0.45 

ft 


ft 
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TABLE  6.  OBSERVED  VERSUS  PREDICTED  WASHOUT  CONCENTRATIONS  - DEA 


Rain 

Collector 

Station 

Measured  Predicted 

Concentration  Concentration 

Cobs.  CH'-  11 

Rain 

Collector 

Station 

Measured  Predicted 

Concentration  Concentration 

^obs.  Si'"  11 

A 8 

0.17  ppm 

A 12 

0.28  ppm 

9 

0.47 

— 

13 

0.41 

— 

20 

0.37 

0.01  ppm 

14 

0.47 

— 

11 

0.43 

0.24 

15 

0.47 

0.001  ppm 

12 

0.57 

0.70 

16 

0.47 

0.019 

13 

0.32 

0.24 

17 

— 

0.112 

14 

0.48 

0.02 

18 

0.63 

0.207 

15 

0.17 

19 

0.56 

0.112 

20 

0.43 

0.019 

21 

0.43 

0.001 

22 

0.42 

C 


high-concentration  data  which  was  extrapolated  over  long  rangea  to  the 
lower  concentrations  and  temperatures  experienced  in  this  study.  A number 
of  factors  could  possibly  render  such  estimates  invalid,  including  the 
effects  of  dissociation,  asaoclatlon,  and  sequestering  by  trace  impurities. 
The  following  sections  of  this  report  describe  the  measurements  performed 
during  the  course  of  this  investigation,  and  the  application  of  these 
results  in  conjunction  with  the  EPAEC  scavenging  model  to  obtain  refined 
washout  estimates. 


CHAPTER  II 


APPARATUS  AND  TECHNIQUE 

Although  a number  of  conventional  techniques  are  available  for  measurement 
of  the  solubility  of  gases  In  liquids  (cf.  HILDEBRAND  and  SCOTT  (1950)), 
the  low-concentration  conditions  of  the  present  study  preclude  the  use  of 
most  such  methods.  The  basic  requirements  of  this  research  are  ones  of 
accurately  setting  and  maintaining  a known  low  concentration  of  solute  In 
either  the  gaseous  (or  aqueous)  phase,  while  accurately  measuring  the 
corresponding  equilibrium  concentration  in  its  aqueous  (or  gaseous) 
counterpart. 

The  basic  approach  of  setting  a known  gaseous-phase  concentration  and 
measuring  Its  aqueous-phase  counterpart  was  attempted  during  the  Initial 
periods  of  this  study.  The  Initial  method  of  maintaining  gaseous-phase 
concentrations  was  the  diffusion-capillary  technique,  as  shown  In  Figure  1. 
With  this  technique  diluent  air  was  passed  at  a known  flow  rate  past  a 
vertical  glass  capillary  containing  a small  amount  of  the  pure  organic,  which 
was  allowed  to  diffuse  Into  the  diluent  stream.  The  dilute  mixture  then 
flowed  Into  a solubility  cell  where  It  was  bubbled  through  a water  sample. 
Upon  approaching  saturation  with  respect  to  the  concentration  of  the  solute 
In  the  gas  stream,  the  sample  was  removed,  and  analyzed  chemically  to  deter- 
mine the  equilibrium  concentration,  and  thus  the  solubility. 

In  theory  this  technique  should  be  totally  adequate  for  the  purpose  of  the 
present  study.  A straightforward  Integration  of  the  duf fusion  equation 
shows  that  the  rate  of  diffusion  of  solute  from  the  glass  capillary  should 
be  given  by 

DC  A 

Flow  rate  (moles /t)  ■ ^ — (”) 

l 

where  D Is  the  diffusion  coefficient  of  solute  In  diluent,  C Is  Its  equl- 

eq 

librium  gas-phase  concentration  (a  function  of  temperature) , A is  the  bore 
area  of  the  capillary,  and  l Is  the  distance  from  the  top  of  the  capillary 


6 


Simplified  Schematic  of  Initial  Solubility  - Measurement  Apparatus 
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to  the  surface  of  the  liquid.  From  equation  (2)  It  is  evident  that  concen- 
tration of  solute  In  the  gaseous  phase  should  be  controllable  through  mani- 
pulation of  a number  of  variables;  and,  provided  that  temperature  and  l 
can  be  maintained  constant,  this  type  of  unit  should  provide  the  long-term 
stabilized  gas  concentrations  required  for  this  study. 

Initial  tests  with  this  apparatus,  however,  proved  unsatisfactory.  It  was 
extremely  difficult  to  obtain  a reliable  calibration  of  the  diffusion  capillary, 
and  highly  erratic  gaseous  concentrations  were  observed  during  solubility 
tests.  After  substantial  amounts  of  testing  it  became  evident  that  pressure 
fluctuations  caused  by  formation  of  bubbles  in  the  solubility  cell  were 
propagating  upstream,  causing  spurious  mixing  In  the  diffusion  capillary. 

Since  there  was  no  straightforward  or  reliable  way  of  countering  this  effect, 
the  diffusion-capillary  approach  was  abandoned  from  further  consideration 
during  this  study. 

In  view  of  this  finding,  a second  dilution  technique,  involving  the  use  of 
permeation  tubes,  was  tested.  Small  quantities  of  the  solute  compounds 
were  encased  lc  polyethylene  vials,  and  located  In  the  diluent  gas  stream 
In  place  of  the  diffusion  capillary.  While  these  units  provided  marginal 
success  in  dampening  short-term  concentration  fluctuations,  they  proved 
unsuccessful  In  providing  the  long-term  stability  required  for  these 
experiments. 

An  additional  difficulty  Involved  with  the  initial  experiments  was  the  long 
time  period  necessary  for  the  approach  to  equilibrium  conditions.  This 
feature  can  be  illustrated  by  rewriting  equation  (1)  in  terms  of  a 
partition  coefficient,  K,  which  relates  gas-  and  liquid-phase  concentrations: 


C,'KC« 

- p/RT 


assuming  ideal-gas  behavior 


(1A) 


Here  C (moles/ml)  Is  the  gas-phase  concentration  of  solute,  rather  than  its 

o 

partial  pressure,  which  appeared  in  (1).  Since 


€ 


8 


(3) 


K - — - — — 

RT  25000 

It  is  obvious  that  extremely  large  volumes  of  gas  must  be  passed  through 
Initially  pure  water  In  the  solubility  cell  before  equilibrium  conditions 
are  approached,  thus  requiring  intolerably  long  experiment  times. 

Because  of  these  difficulties  a reverse  experiment,  based  upon  fixing  known 
concentrations  of  solute  in  the  aqueous  phase  and  measuring  equilibrium 
gas-phase  concentrations  was  tested.  Shown  In  Figure  2,  this  experiment 
involved  bubbling  ultrapure  helium  (supplied  by  the  U.  S.  Bureau  of  Mines) 
from  a regulated  cylinder  through  a flow  meter,  a heat-transfer  coll,  and 
finally  a solubility  cell,  containing  a known  solution  of  solute  In  triple- 
distilled,  organic-free  water.  Helium,  thus  saturated  with  solute,  was 
conducted  into  a gas-chromatograph  sampling  valve  where  It  was  subsequently 
Injected  into  a chromatograph  for  analysis.  The  chromatograph  was  a flame 
Ionization  unit,  manufactured  by  Analytical  Instruments  Development,  Inc. 

Owing  to  the  high  purity  of  both  the  solvent  gas  and  the  water,  the  solute 
vapor  was  the  ouly  organic  In  the  chromatograph  sample,  and  thus  the  separa- 
tion requirements  were  minimal.  All  that  was  required  was  separation  of 
the  solute  peak  from  the  minimal  air  peak,  and  this  was  accomplished  using 
an  unpacked  1/8"  stainless  steel  column,  approximately  100  cm  long.  A 
typical  peak  from  an  ethyl  acetoacetate  sample  la  shown  in  Figure  3. 

In  contrast  to  the  previous  experiments,  the  behavior  discussed  In  conjunc- 
tion with  equations  (1A)  and  (3)  is  one  sense  advantageous  for  the  present 
technique.  Since  the  partition  coefficient  la  small  for  both  ethyl  aceto- 
acetate and  diethylamine,  comparitively  large  quantities  of  divent  gas  can 
be  passed  through  the  solubility  cell  with  negligible  change  in  liquid 
phase  concentration.  This  results  In  a relatively  stable  system,  \ rovidlng 
ample  time  for  decay  of  transients  prior  to  final  measurement. 
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Figure  2.  Schematic  of  Final  Solubility  — Measurement  Apparatus 


In  another  sense,  however,  this  behavior  is  disadvantageous.  Gas-phase 
concentrations,  rather  than  those  in  the  aqueous  phase,  are  measured  in  this 
approach;  and  because  of  the  low  value  of  K this  required  the  observation  of 
comparitively  low  concentrations,  which  are  more  difficult  to  measure.  Owing 
to  this  feature  solubilities  at  gas-phase  concentrations  below  a few  tenths 
of  a part  per  million  could  not  be  determined  using  this  technique. 

Additional  aspects  of  this  experimental  method  give  rise  to  other  potential 
measurement  errors.  The  method  depends  upon  the  assumption  that  the  gas 
stream  leaving  the  liquid  is  saturated  with  solvent,  and  this  in  turn 
requires  that  the  carrier  gas  be  passed  through  the  system  at  a rate 
sufficiently  slow  to  ensure  that  mass  transfer  from  the  liquid  to  the  gas 
had  proceeded  essentially  to  completion.  In  addition,  too  high  a gas  flow 
rate  could  possibly  entrain  liquid  droplets,  and  deposit  them  in  the  line 
between  the  solubility  cell  and  the  chromatograph  sampling  valve,  where 
they  might  desorb  additional  solute  and  thus  produce  high  results. 

Both  of  the  above  effects  were  examined  by  varying  the  gas  flow  rate  through 
a test  mixture  and  observing  the  resulting  concentrations  of  solute  in  the 
gas  stream.  From  these  tests  it  was  determined  that  gas  concentration  was 
virtually  independent  of  flow  rate  for  the  solubility  cell  employed  (250  ml 
Smith-Greenburg  impinger  containing  40  m.  liquid),  at  flow  rates  below  about 
40  ml/minute.  Above  40  ml/minute  higher  concentrations  were  observed, 
ostensibly  from  entrainment  of  liquid  droplets.  For  this  reason,  all 
solubility  experiments  were  conducted  at  nominal  flow  rates  of  10  ml/minute. 

A typical  solubility  experiment  was  conducted  as  follows:  A solution  of 

known  solvent  concentration  was  prepared  by  volumetric  dilution  in  triple- 
distilled  water,  and  40  ml  of  this  solution  was  placed  in  the  solubility 
cell.  The  cell  was  in  turn  placed  in  the  thermostatic  bath  (normally  at 
5°C)  and  allowed  to  approach  thermal  equilibrium.  Upon  approaching  the 
bath  temperature  the  gas  flow  was  turned  on,  and  the  resulting  concentra- 
tion of  solute  in  the  gas  was  monitored  until  all  transients  had  passed, 
and  the  true  equilibrium  concentration  was  attained. 
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The  gas  chromatograph  was  calibrated  each  day  that  solubility  measurements 
were  performed.  This  was  accomplished  using  a syringe  pump,  which  injected 
an  air-solute  mixture  into  a constant  flow  of  helium  at  selected  flow 
rates.  Known  concentrations  of  solute  in  air  were  loaded  into  the  syringe 
by  withdrawing  volumes  of  vapor  at  equilibrium  with  the  liquid  phase  at  a known 
temperature  from  an  expanding- volume  vessel.  Vapor  pressures  of  the 
solutes  were  calculated  for  this  purpose  from  the  expressions 


—6098 

Jin  p * — m — + 20.21  for  ethyl  acetoacetate  » (4) 

and 

£n  p =■  + 18.089  for  diethylamine  , (5) 


where  p is  vapor  pressure  in  mm  Hg,  and  T is  in  degrees  Kelvin.  Plots  of 
the  calibration  curves  produced  in  this  manner,  along  with  the  numerical 
data,  are  given  in  the  Appendix. 


The  flow  rate  of  diluent  gas  was  monitored  with  a rotameter,  which  had 
been  calibrated  using  a wet  test  meter.  Flow  rates  from  the  syringe 
pump  were  calibrated  using  a soap-film  flowmeter.  All  tubing  downstream 
from  the  mixing  points  of  both  the  calibration  apparatus  and  solubility 
cell  was  1/A  inch  Teflon,  to  minimize  sorption  of  the  organic  vapors. 
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The  experimental  results  of  this  study  are  given  in  Tables  7 and  8.  These 
Include  the  data  obtained  for  both  solutes  at  5.0°C,  in  addition  to  data 
for  higher  temperatures.  Data  at  5.0#C  for  acidified  water  solutions  (HC1 
added)  for  diethylamine  are  given  as  well.  These  results  will  be  discussed 
in  the  following  chapter. 


TABLE  7.  SUMMARY  OF  SOLUBILITY  MEASUREMENTS  FOR  DIETHYLAMINE 


Solution 

Molarity 

Solute  Concentration 

Solute  Mixing 

Temperature 

Added 

Acid 

in  lluqid 

Ratio  in  Gas 

®K 

molar 

moles/ml 

pptn 

278.2 

0 

1.22 

X 

Iff6 

2.7 

278.2 

0 

2.43 

X 

l<f6 

6.0 

278.2 

0 

4.87 

X 

o 

1 

23.7 

278.2 

0 

• 

OO 

X 

JO'6 

24.2 

278.2 

0 

9.73 

X 

io"5 

41.4 

278.2 

0 

9.63 

X 

ID’5 

355 

278.2 

2.50 

X 

*-• 

o 

i 

u> 

2.37 

X 

icf6 

1.5 

278.2 

2.50 

X 

10"  3 

4.75 

X 

VO 

1 

o 

H 

6.4 

278.2 

2.50 

X 

l(f3 

9.49 

X 

vO 

1 

O 

H 

32.7 

278.2 

2.50 

X 

10-3 

9.40 

X 

IO-5 

347 

279.9 

0 

9.73 

X 

*-* 

o 

1 

CT* 

43.2 

282.1 

0 

9.73 

•x 

io-6 

44.8 

291.6 

0 

9. 73 

X 

Hf6 

137 

o 
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TABLE  8.  SUMMARY  OF 

SOLUBILITY  MEASUREMENTS 

FOR  ETHYL  ACETOACETATE 

Solution 

Solute  Concentration 

Solute  Mixing 

£ 

Temperature 

In  Liquid 

Ratio  in  Gas 

•k 

moles /ml 

ppm 

278.2 

7.85  x 10"6 

4.7 

c 

278.2 

3.94  x 10“6 

3.1 

278.2 

7.88  x 10“5 

12.7 

278.2 

7.88  x 10"5 

12.3 

278.2 

7.16  x 10“6 

3.2 

0 

278.2 

3.75  x 10“6 

2.7 

278.2 

1.92  x 10"6 

2.0 

269.2 

7.16  x 10"6 

11.0 

Q 
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CHAPTER  IV 


DISCUSSION  OF  RESULTS 


EXPERIMENTAL  ERRORS 

Accuracy  of  the  experimental  results  was  principally  dependent  upon  the 
accuracy  of  the  gas-phase  concentration  measurements.  These  measurements 
depended,  In  turn,  upon  the  calibration  of  the  gas-chromatograph,  and  the 
degree  of  error  Involved  In  reading  chromatograph  output.  The  minimum 
detectable  gas-phase  concentrations  were  of  the  order  of  a few  tenths  of  a 
part  per  million  for  both  diethylamlne  and  ethyl  acetoacetate,  with  ethyl 
acetoacetate  exhibiting  somewhat  higher  sensitivity.  At  higher  concen- 
trations, serial  measurements  Indicated  a reporduclblllty  of  about  ± 5Z. 
These  features  give  rise  to  an  overall  repeatability*  of  the  measurements 
corresponding  to  the  equations 


and 


E ■ 0.5  + .05  x PP®  for  diethylamlne 
E - 0.3  + .05  x PP®  for  ethyl  acetoacetate 


(6) 

(7) 


where  E Is  the  expected  error,  and  x is  the  measured  mixing  ratio  In  parts 
per  million. 


Calibration  of  the  gas  chromatograph  Involved  the  same  peak-reading  errors 
as  the  sampling  procedure.  Serial  analyses  of  calibration  mixtures 
Indicated  a similar  repeatability;  thus  a conservatively  large  error 
estimate  Is  obtained  for  the  combination  of  the  effects  simply  by  doubling 
the  above  expressions 


Etqt  1 +0.1  x PP®  for  diethylamlne  (8) 

and 

ET0T  " + X PP®  f°r  ethyl  acetoacetate  . (9) 

As  stated  previously,  additional  errors  arising  from  temperature  control, 
flow  control,  standards  preparation,  and  entrainment  and  disequilibrium 

*using  a consistent  calibration  standard. 
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effects  were  comparltlvely  minor;  hence  (8)  and  (9)  should  represent  the 
total  expected  error  reasonably  well  except  for  the  case  of  low-concentra- 
tion ethyl  acetoacetate,  which  was  found  to  vary  radically  with  concentra- 
tion and  possibly  lead  to  higher  errors.  This  aspect  will  be  discussed 
later  In  this  chapter,  where  this  compound  is  addressed  more  specifically. 


DIETHYLAMINE 


In  analyzing  the  diethyiamlne  solubility  data  it  Is  helpful  to  set  forth 
a mechanism  for  Interaction  of  this  material  with  water.  Ostensibly,  this 
Interaction  consists  of  an  absorption  from  the  gas  phase  with  a subsequent 
Ionization  In  water  solution.  Representing  these  steps  In  the  conventional 
manner,  we  have 


(c2h5)2  NH  | 

H20  + (C2H5)2  NH 
where  the  equilibrium 

Ko 

*b 


K 


gas  t'  (C2H5^2  1,11  I water 

+ 

| water  -p2— -»  (C2H5>2  NH2  + 0H“ 

constants  Kq  and  are  given  by 

I(C2H5)2  NH  | water] 

" [(C2H5)2  NH  | gas] 

[(C2H5)2  NH2]  [0H_] 

“ [(C2H5)2  NH  | water]  * 


(10) 

(11) 


(12) 


(13) 


Since  total  dissolved  diethyiamlne  exists  as  two  species,  l.e., 

DE^  - [(C2H5)2  N2H+]  + l(C2H5)2  NH  | water]  (14) 

one  can  formulate  an  expression  relating  gaseous  to  dissolved  solute 

DEA_  [OH-] 

[(C2H5)2  NH  | gas]  - Kq  (1^  + [OH"])  * 

An  accompanying  expression  for  [OH-]  can  be  obtained  by  employing  the 
'equilibrium  relationship  for  dissociation  of  water 
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with  the  Ion  balance  for  the  solution 


'“"total1  ■ t0HV  + '""'dea1  • (l7) 

The  result  Is 

IOHTOTAlJ  + (H3°  don  + ^OHTOTAL^  + 

(18) 

(-Kw  DEAt  + H30+don  V [0H-total]  - K,  Kw  - 0 . 

which  expresses  the  OH  concentration  from  all  sources,  f0HxoTAL^’  as  a 
cubic  equation.  In  the  above  equations,  OH  w and  OH  D£A  denote  OH  arising 
from  Ionization  of  water  and  diethylamine,  respectively.  H3°+,jon  is  the 
hydrogen  ion  contributed  to  the  system  from  sources  other  than  dissociation 
of  water;  this  term  accounts  for  the  influence  of  acid  on  the  solubility  of 
diethylamine . 

If  values  of  the  equilibrium  constants  were  known,  one  could  proceed 

Immediately  to  solve  (15)  and  (18)  to  obtain  a family  of  curves  relating 

the  partial  pressure  of  diethylamine  above  water  solutions.  Indeed,  values 

for  K are  well  established  (RIEMAN,  NEUSS,  and  NAIMAN  (1951))  and  a few 
w 

measurements  of  are  available  in  the  literature,  which  can  be  extrapolated 
to  temperature  conditions  of  interest  with  good  reliability.  Utilizing  the 
data  of  Evans  and  Hayamann  (1951)  for  this  purpose  in  conjunction  with  the 
van't  Hoff  expression,  one  obtains  the  form 

fcn  ^ - - 8.4726  + 21-1-^-15  . (19) 

Solving  equation  (18)  for  values  of  Kw  and  appropriate  to  5®C  (1.86  x 10-^ 
and  4.48  x 10  *,  respectively)  the  results  given  in  Table  9 are  obtained. 


TABLE  9.  ROOTS  OF  EQUATION  18  FOR  VARIOUS  CONDITIONS 


Dissolved  Diethylamine 

OH  Concentration 

OH  Concentration 

Concentration 

H,0+,  - 0 

3 don 

H.O+  - 

3 don 

2.5  x 10“ 

molar 

molar 

molar 

1 X 

10"3 

4.82  x 10"A 

1.24  x 

i<r12 

5 x 

10"  3 

1.29  x 10“ 3 

3.41  x 

i<f4 

1 X 

10"2 

1.90  x 10"3 

8.78  x 

i<f4 

5 x 

10"2 

4.51  x 10"3 

3.37  x 

lo-3 

1 X 

10"1 

6.47  x 10~3 

5.30  x 

10-3 

Since  Kq  Is  the  only  remaining  parameter  In  equations  (15)  and  (18)  it  can 
be  adjusted  to  fit  the  experimental  results.  If  the  mechanism  set  forth  in 
this  section  is  correct,  then  the  experimental  data  should  conform  to  the 
solutions  to  (15)  and  (18)  over  the  complete  range  of  variables. 

A visual  fit  of  the  data  to  these  equations  suggested  that  Kq  at  5°C 
should  be  approximately  equal  to  5060,  and  a corresponding  comparison  of 
the  experimental  and  theoretical  results  is  given  in  Figure  4.  Here  it 
should  be  emphasized  again  that  the  curves  were  adjusted  to  the  data  by 
choosing  an  appropriate  value  of  Kq.  The  good  fit  to  all  of  the  data 
points  over  the  ranges  of  concentration  and  acidity  investigated,  however, 
suggests  that  the  postulated  mechanism  is  a valid  one,  and  provides  a 
basis  for  the  calculation  of  solubility  under  all  low-concentration 
conditions. 

ETHYL  ACETOACETATE 

In  contrast  to  the  basic  behavior  of  diethylamine,  ethyl  acetoacetate 
exhibits  mildly  acedlc  behavior  in  water  solution.  Its  small  dissociation 
constant,  however,  of  the  order  of  2 x 10  ^ (EIDINOFF  (1945)),  limits  the 
significance  of  dissociation  phenomena  insofar  as  solubility  is  concerned, 
except  possibly  at  extremely  low  concentrations  or  under  highly  basic 
conditions.  This  compound  also  enters  into  a tautomer  it  at  ion  reaction  and 


exists  in  both  keto  and  enol  forms,  which  have  distinctly  different  chemical 
characteristics  (cf.  MORRISON  and  BOYD  (1959),  RIDDICK  and  BUNGER  (1970)). 

At  room  temperature  ethyl  acetoacetate  contains  about  eight  percent  of  the 
enol  form. 

Because  of  the  noted  complexity  of  the  ethyl  acetoacetate  system,  it  is  most 
expedient  here  to  express  the  solubility  data  presented  in  the  previous 
chapter  directly  in  terms  of  the  solubility  parameter  H' . This  is  done  in 
the  graph  in  Figure  5 where  H'  is  plotted  as  a function  of  solute  concen- 
tration in  the  liquid  phase.  Of  particular  interest  here  is  the  rapid 
decrease  of  solubility  (increase  of  H')  at  concentrations  lower  than  about 

-3 

10  molar.  This  feature  made  accurate  measurements  of  solubility  more 
difficult  under  such  conditions,  since  slight  changes  of  solute  concentra- 
tion had  pronounced  effects  on  the  corresponding  values  of  H' . 

It  is  difficult  to  explain  the  phenomena  exhibited  in  Figure  5.  Although 
the  complex  chemistry  of  ethyl  acetoacetate  may  be  the  cause  of  this 
behavior,  it  seems  more  likely  that  it  is  based  in  the  strongly  self- 
associated  nature  of  this  compound.  This  type  of  behavior  is  conventionally 
associated  with  positive  deviations  from  Raoult's  Lav  (DENBIGH  (1957)),  and 
rests  simply  in  the  fact  that  clusters  of  solute  molecules  can  exist  at 
considerably  lower  energy  than  solvent-solute  associated  molecules.  At 
very  low  concentrations,  where  solute-solute  interactions  become  less 
likely,  the  resulting  high  energy  state  will  be  reflected  in  an  Increase 
in  solute  vapor  pressure. 

Such  behavior  implies  that  a limiting  value  of  H'  will  be  approached  at 
concentrations  approaching  zero.  This  value  was  determined  from  the  5*C 
data  by  plotting  the  logarithm  of  H'  versus  concentration  and  linearly 
extrapolating  to  zero: 


H'  ■ 1.4  atm.  ml/mole  . 
o 

Combining  this  with  a visual  fit  to  the  low- concent  rat  ion  data  provides 
the  empirical  form  for  the  vater  solubility  of  ethyl  acetoacetate  at  5*C: 
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Figure  5.  Experimental  Solubilities  for  Ethyl 


(20) 


H'  ■ 1.4  - 1.68  x 10“  C atm.  ml/mole 
(0  < C < 5 x 10"6) 

where  C Is  solute  concentration  In  the  liquid  phase  In  moles/ml. 

TEMPERATURE  DEPENDENCE 

Although  the  bulk  of  data  obtained  In  this  project  were  acquired  at  5*C. 
additional  results,  obtained  at  higher  temperatures,  allow  estimates  of 
temperature  dependence  of  the  solubilities  of  both  diethylamine  and 
ethyl  acetoacetate.  For  diethylamine,  the  temperature  dependence  of  both 
and  are  known;  thus  if  that  of  Kq  can  be  determined  total  temperature 
dependence  can  be  assessed. 

Applying  (IS)  and  (18)  to  the  experimental  data  point  for  291.6  *K  In 
Table  7 results  in  a corresponding  value  for  Kq  of  approximately  1380. 
Utilizing  an  equation  of  the  form  (19)  and  fitting  the  278.2  *K  and 
291.6  *K  results  gives 

In  Kq  • - 19.75  + . (21) 

A similar  treatment  can  be  used  to  derive  an  analagous  equation  for  ethyl 
acetoacetate,  although  this  must  be  considered  highly  approximate  In  view 
of  the  complicated  behavior  of  this  substance  discussed  previously.  Using 
the  data  points  at  7.16  x 10  ^ moles/cc  at  278.2  and  296.2  degrees  K (H* 
of  .45  and  1.54  atm.  ml/mole,  respectively)  In  conjunction  with  an  equation 
of  the  form  (19)  gives 

An  H'  - 19.446  - —2.  (22) 

at  7.16  x 10  ® moles/cc.  Changes  In  H'  at  other  concentrations  can  be 
estimated  to  a first  approximation,  at  least,  by  assuming  proportionate 
changes  with  changes  In  temperature. 
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REVISED  SCAVENGING  CALCULATIONS 

The  measurements  obtained  In  this  study  showed  the  solubilities  of  diethyl- 
amine  and  ethyl  acetoacetate,  respectively,  to  be  higher  and  lower  than 
estimated  from  previous  measurements.  As  mentioned  in  the  first  chapter, 
this  is  consistant  with  the  comparison  of  predicted  and  measured  washout 
rates  of  these  compounds,  and  it  is  of  interest  here  to  perform  revised 
estimates  based  upon  the  newly-acquired  solubility  data. 

Such  estimates,  based  upon  EPAEC  calculations  in  conjunction  with  the 

solubility  equations  (15).  (18),  and  (20)  are  shown  in  Tables  10-14,  which 

may  be  compared  directly  with  Tables  2-6  in  Chapter  I.  From  these  results 

the  following  observations  can  be  made: 

o estimates  based  on  the  revised  solubility  data  are  in 
significantly  better  agreement  with  observation  than 
previous  estimates, 

o better  agreement  between  observation  and  theory  occurs  for  the  more 

distant  sampling  line  (line  B)  than  for  the  inner  sampling  line  (line  A) , 

o Dlethylamlne  measurements  exhibit  closer  agreement  with 
theory  than  do  those  for  ethyl  acetoacetate. 

In  the  case  of  ethyl  acetoacetate  it  is  readily  apparent  that,  although  the 
modification  of  the  solubility  estimates  did  enhance  the  agreement  between 
observation  and  theory,  calculated  concentrations  typically  exceed  experi- 
mental values  by  a factor  of  five  or  more.  This  behavior  indicates  that 
scavenging  is  Influenced  by  additional  rate  phenomena,  and  suggests  that 
liquid-phase  mixing  within  the  raindrops  may  be  emerging  as  an  Important 
factor  in  this  respect. 


Liquid-phase  mixing  is  dependent  upon  molecular  and  convective  transport 
processes  within  falling  raindrops,  and  its  influence  should  increase  as  the 
diffusion  coefficient  and  solubility  of  the  material  become  small.  Since  no 
dlffuslvity  data  are  available  for  the  ethyl  acetoacetate-water  system  this 
effect  is  difficult  to  assess.  Assuming  a moderately  low  diffusion  coeffic- 
ient for  this  material,  however,  it  is  not  unreasonable  to  expect  a factor 
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TABLE  10.  REVISED  PREDICTIONS  OF  WASHOUT  CONCENTRATIONS 


* 


e 


Experiment  1 (EAA) 


Z 


c 


c 


o 


o 


c 


o 


Rain 

Collector 

Station 

Measured 

Concentration 

C . 
obs. 

Predicted 

Concentration 

C 

Rain  Measured 

Collector  Concentration 

Station  C . 

obs. 

Predicted 

Concentration 

C 

A 3 

<0.10  ppm 

B 8 

— 

— 

4 

0.16 

0.03  ppm 

9 

— 

0.81  ppm 

5 

0.50 

0.27 

10 

— 

1.48 

6 

1.1 

1.55 

11 

1.0  ppm 

2.15 

7 

1.5 

6.34 

12 

0.70 

2.46 

8 

1.6 

14.43 

13 

1.2 

2.98 

9 

1.5 

6.34 

14 

0.78 

2.46 

10 

1.2 

1.55 

15 

0.73 

2.15 

11 

1.0 

0.27 

16 

0.33 

1.48 

12 

0.32 

0.03 

1' 

0.35 

0.81 

13 

0.1 

— 

18 

0.12 

— 

TABLE  11.  REVISED  PREDICTIONS  OF  WASHOUT  CONCENTRATIONS 
Experiment  2 (EAA) 

Rain  Measured  Predicted  Rain  Measured  Predicted 

Collector  Concentration  Concentration  Collector  Concentration  Concentration 
Station  Cob>_  c Station  c 

A 1 

<•10  ppm 

0.22  ppm 

B 5 

<0.10  ppm 

0.69  ppm 

2 

<•10 

0.69 

6 

0.26 

1.0 

3 

0.80 

1.8 

7 

0.16 

1.3 

4 

1.0 

3.6 

8 

0.31 

1.6 

5 

1.5 

5.9 

9 

0.36 

1.8 

6 

2.0 

9.4 

10 

0.31 

1.9 

7 

1.6 

5.9 

11 

0.27 

2.0 

8 

0.77 

3.6 

12 

0.37 

1.9 

9 

0.20 

1.8 

13 

0.28 

1.8 

10 

<0.10 

0.69 

14 

0.12 

1.6 

11 

<0.10 

0.22 

15 

<0.10 

1.3 

12 

<0.10 

0.06 

16 

<0.10 

1.0 

13 

<0.10 

17 

<0.10 

0.69 

18 

<0.10 

c 
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TABLE  12.  REVISED  PREDICTIONS  OF  WASHOUT  CONCENTRATIONS 


Experiment  3 (DEA) 


Rain 

Measured 

Predicted 

Rain 

Measured 

Predicted 

Collector 

Concentration 

Concentration 

Collector 

Concentration 

Concentration 

Station 

C . 
obs. 

C 

Station 

^obs. 

C 

A 4 

0.6  ppm 

«... 

B 11 

0.9  ppm 

0.7  ppm 

5 

0.7 

— 

12 

1.7 

1.4 

6 

2.1 

— 

13 

2.2 

2.0 

7 

1.4 

0.2  ppm 

14 

2.7 

2.5 

8 

1.3 

1.4 

15 

3.0 

2.8 

9 

2.1 

4.9 

16 

2.1 

2.5 

10 

2.5 

7.3 

17 

2.8 

2.0 

11 

1.0 

4.9 

18 

1.1 

1.4 

12 

0.6 

1.4 

19 

1.2 

0.7 

13 

<0.1 

0.2 

14 

<0.1 

— 

TABLE  13.  REVISED  PREDICTIONS 

1 OF  WASHOUT  CONCENTRATIONS 

Experiment 

4 (DEA) 

Rain 

Measured 

Predicted 

Rain 

Measured 

Predicted 

Collector 

Concentration 

Concentration 

Collector 

Concentration 

Concentration 

Station 

C . 
obs. 

C 

Station 

C . 
obs. 

C 

A 6 

.19  ppm 

.1  ppm 

B 12 

.21  ppm 

— - 

7 

.55 

.01 

13 

.42 

.15  pom 

8 

.62 

.46 

14 

1.3 

1.0 

9 

.96 

5.2 

15 

1.3 

2.1 

10 

1.2 

11.2 

16 

1.7 

3.3 

11 

1.0 

5.2 

17 

2.2 

4.6 

12 

.85 

.46 

18 

1.7 

3.3 

13 

.56 

.01 

19 

1.4 

2.1 

14 

.39 

— 

20 

1.8 

1.0 

21 

1.3 

.15 

22 

.80 

— 

23 

.45 

O 
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TABLE  14.  REVISED  PREDICTIONS  OF  WASHOUT  CONCENTRATIONS 


Experiment  5 (DEA) 


Rain 

Measured 

Predicted 

Rain 

Measured 

Predicted 

Collector 

Concentration 

Concentration 

Collector 

Concentration 

Concentration 

Station 

C . 
obs. 

C 

Station 

C . 

obs. 

C 

A 8 

0.17  ppm 

B 12 

0.28  ppm 

.1  ppm 

9 

0.47 

.02  ppm 

13 

0.41 

.2 

10 

0.37 

.41 

14 

0.47 

.8 

11 

0.43 

3.1 

15 

0.47 

1.3 

12 

0.57 

6.0 

16 

0.47 

1.7 

13 

0.32 

3.1 

17 

— 

2.2 

14 

0.48 

.41 

18 

0.63 

1.7 

15 

0.17 

.20 

19 

0.56 

1.3 

20 

0.43 

.8 

21 

0.43 

.2 

22 

0.42 

.1 

If 
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of  ten  reduction  In  washout  concentrations  from  this  effect.  Quantitative 
assessments  of  the  effects  of  liquid  mixing,  if  diffusion  coefficients  are 
known,  can  be  performed  by  the  EPAEC  model  (DANA,  HALES,  and  WOLF  (1973)). 

Scavenging  calculations  for  dlethylamlne  are  complicated  by  the  presence  of 
atmospheric  COj,  which  forms  acid  upon  dissolution  in  water  and  thus  modifies 
dlethylamlne  solubility.  Conversely,  dissolved  dlethylamlne  modifies  the 
solubility  of  COg  In  water,  and  the  simultaneous  presence  of  these  materials 
should  enhance  the  scavenging  rates  of  both.  The  EPAEC  model  cannot  account 
for  scavenging  of  multiple,  reactive  components*,  and  for  this  reason  we  have 
chosen  to  account  approximately  for  the  Influence  of  CO,  by  choosing  a nominal 
constant,  associated  value  of  [H^O  ^on]  - 3.2  x 10  molar.  This  approximation 
should  be  reasonably  valid  for  this  purpose  owing  to  the  relatively  slow 
dissociation  rate  of  CO^  in  wafer,  and  to  the  rather  secondary  effects  of 
hydrogen  Ion  on  solubility  In  this  concentration  range. 

As  seen  from  the  data  In  Tables  18-20,  the  computed  washout  concentrations 
for  the  innermost  arc  appear  to  be  somewhat  higher  and  more  constricted  In 
the  crosswind  direction  than  those  actually  observed.  Such  behavior  Is 
typical  of  similar  measurements  for  SO2  under  such  conditions,  and  Is  thought  to 
be  caused  primarily  by  undercutting  by  the  rain  and  meandering  by  the  true 
plume,  which  Is  approximated  as  a point  source  In  the  EPAEC  model.  The 
agreement  between  experiment  and  theory  on  the  outer  sampling  line  Is 
excellent;  this  Is  undoubtedly  caused  in  part  by  fortuoslty,  but  Is  considered 
also  to  be  strong  supporting  evidence  for  this  modeling  approach.  From  these 
Investigations  It  Is  concluded  that  although  a number  of  microphysical  phenomena 
may  emerge  as  rate-influencing  steps  In  the  gas-scavenging  process,  the 
modeling  approach  presented  in  this  and  the  previous  report  provides  a 
reasonable  and  convenient  means  for  generalized  gas-scavenging  analysis. 


*/  new  code.  Scavenging  Model  Incorporating  Chemical  Kinetics  (SMICK)  has 
recently  been  developed  at  this  laboratory  for  this  purpose  under  contract 
to  the  EPA. 
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APPENDIX  A 

CHROMATOGRAPH  CALIBRATIONS  AND  RAW  DATA 


DIETHYLAMINE 


PPM  PEA 
4.43 
8.87 

17.7 
44.3 

88.7 
177 


Chromatograph  Calibrations 
9 July  75 

PEAK  HT 
29 
85 
240 
1190 
3170 
9290 


225 

450 


11  July  75 


9980 

21800 


Solubility  Measurements 


Liquid  Cone, 
mole 8 /ml 


1.22  x 10 

7.43  x 10"6 

2.37  x 10"6 

4.87  x 10"6 

4.75  x 10”6 

4.87  x 10’6 

9.73  x 10"6 

9.49  x 10'6 

9.73  x 10"6  (6,7*C) 

9.73  x 10"6  (8,9*0 

9.73  x 10'6  (18,4*0 


'"don1 

molar 

PEAK 

0 

12 

0 

43 

2.5  x 10"3 

5 

0 

378 

2.5  x 10“3 

48 

0 

392 

0 

917 

2.5  x 10“3 

632 

0 

982 

0 

1040 

0 

6080 
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Solubility  Measurements  (continued) 


Liquid  Cone. 

[sd<J 

PEAK  H‘ 

moles/ml 

molar 

9.63  x 10"5 

0 

17200 

9.40  x 10-5 

2.5  x *0-3 

16900 

ETHYL  ACETOACETATE 


Chromatograph  Calibrations 
6 July  75 


PPMEAA 


PEAK  HT 


1.95 

3.90 

7.80 

15.6 


55 

136 

288 

708 


11  July  75 


0.74  13 
1.48  38 
2.96  102 
5.92  248 
11.8  610 
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Solubility  Measurements 


O 


0 


c 


0 


c 


Liquid  Cone, 
moles/ml 


6 July  75 

7.88  x 10"6 
3.94  x l(f6 


11  July  75 

7.88  x 10-5 
7.88  x 10'5 
7.16  x l(f6 
3.75  x 10“6 
1.92  x 10"6 
7.16  x l(f6  (23.0°C) 


O 


Q 


G 


C 


t 


32 


PEAK  HT 


162 

98 

720 

690 

112 

88 

60 

592 


